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ABSTRACT 

We discuss the fragmentation of primordial gas clouds in the universe after 
decoupling. Comparing the time scale of collapse with that of fragmentation, 
we obtain the typical mass of a fragment both numerically and analytically. It 
is shown that the estimated mass gives the minimum mass of a fragment which 
is formed from the primordial gas cloud and is essentially determined by the 
Chandrasekhar mass. 



Subject headings: galaxies: formation — molecular processes — stars: formation 



1. Introduction 



The first stars formed after decoupling had no metals heavier than lithium, so that 
star formation proceeded without grains. This means that the isothermal approximation, 
which is often used to investigate star formation in the present molecular cloud (its validity 
is explained by Hayashi & Nakano 1965), is not applicable and that star formation in the 
primordial gas cloud may be completely different. Up to now many authors have studied 
dynamical and thermal evolutions of primordial gas clouds (Matsuda, Sato, & Takeda 1969; 
Yoneyama 1972; Hutchins 1976; Yoshii & Sabano 1979; Carlberg 1981; Palla, Salpeter, 
& Stahler 1983; Lahav 1986; Puy & Signore 1996; Susa, Uehara, & Nishi 1996, hereafter 
paper I). One of the most important quantity in such studies is the minimum Jeans mass in 
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the primordial gas cloud in relation to the initial mass function. In a spherical free-falling 
cloud, however, the density perturbation grows only in a power law (Hunter 1964) as in the 
expanding universe, so that the minimum Jeans mass does not necessarily mean the mass 
of the fragment (Larson 1985 and references therein). In general a collapsing cloud does 
not necessarily fragment at the instant when the Jeans mass reaches the minimum, but 
fragments when the collapse of the cloud is almost halted (see §|3|). 

It is known that a spherical cloud in pressure-free collapse is unstable against 
non-spherical perturbations (Lin, Mestel, & Shu 1965; Hutchins 1976; paper I). This means 
that the collapsing primordial gas cloud is not spherical even without angular momentum 
in general. Therefore the evolution of a non-spherical cloud must be considered. In one 
of the plausible scenarios, shock waves occur in the collapsing non-spherical primordial 
gas cloud. In the post shock flow, due to cooling by hydrogen molecules formed through 
non-equilibrium recombination, a cooled sheet (T <^ 500K) is formed (MacLow & Shull 
1986; Shapiro & Kang 1987, hereafter SK). This cooled sheet fragments more easily into 
filaments (Miyama, Narita, & Hayashi 1987a, 1987b) than into spherical clouds. As the 



virial temperature of a filament is essentially constant (see eq. [ II]), its evolution is much 
different from that of a spherical cloud whose virial temperature increases as the cloud 
contracts. 

In this Letter we investigate the dynamical and thermal evolution of an infinitely 
long cylinder to model a filament formed in the post shock cooled sheet and estimate the 
minimum mass of the fragment of the cylinder. The formulation of our model is shown in 
§0. In §|3] we show how to estimate the mass of the fragment in our model. The numerical 
results are shown in §f|. We analytically estimate the final fragment mass in will be 

devoted to discussions . 



2. Basic equations 
2.1. Equation of motion of the cloud for cylindrical collapse 

For an isothermal cylinder there is a critical line density M C (T) (Ostriker 1964)given by 

M C (T) = (1) 

where T, /z, and me represent the isothermal temperature, the mean molecular weight and 
the mass of hydrogen atoms, respectively. An isothermal cylinder with line density M 
larger than M C (T) has no equilibrium and collapses, while an isothermal cylinder with M 
less than M C (T) expands if the external pressure does not exist. 
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When the cylindrical cloud is formed as a result of the fragmentation of the sheet, the 
gas pressure is comparable to the gravitational force. So the effect of the gas pressure must 
be included in the equation of motion of the cylindrical cloud. To know the dynamical 
evolution of such a cylinder we begin with the virial equation for a cylinder, 

)^± = 2T + 2n-GM\ (2) 

where 

1=1 pr 2 dV, T= f \pu 2 dV, 77 = [ pdV, (3) 
Jv Jv 2 Jv 

and the integration is effected over the volume per unit length of the cylinder, V. Here 

we approximate that the cylinder is uniform and has constant density p, scale radius 



R = ^M/irp, and temperature T. Then above 7, T, and 77 become 

nn ^,JdR\ 2 „ k B T 



MR , T = -M[— , 77 = ^M. (4) 
2 4 \ dry pm^ 



Substituting (|]) in equation (0) we obtain 



dt 2 R 

We describe the evolution of the cylinder by this equation. 



t]1U ^(M-M C (T)). (5) 



2.2. Energy equation and chemical reactions 

The equations describing thermal processes are the same as those of paper I. The time 
evolution of the cloud temperature is described by the energy equation given by 

de d 1 

—T = ~P~r, A-rad — J^chem , (6) 

dt dt p 

where e, A ra d and A c h em are the thermal energy per unit mass , the cooling by the radiation 
and the cooling/heating by chemical reactions, respectively. In A rad we take into account 
H 2 line emissions due to the vibrational and rotational transitions between excited levels. 
The gas pressure P is evaluated by the equation of state of the ideal gas. The details of 
Arad are shown in Appendix B of paper I and A c h era is shown in SK. 

We use the same chemical reaction rates as Palla et al. (1983). We consider 
H, H2, H + , H~, and e~ as chemical compositions of the cloud. Helium and its ions are not 
included in our calculations because they are chemically inert at such low temperatures 
(T 10 4 K) as considered here. Hereafter we use fi to denote the fractional number 
abundance of the i-th composition. 
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3. The epoch of fragmentation and the mass of a fragment 

The dispersion relation of the linear density perturbations of an isothermal sheet with 
surface density X shows that the perturbation with wavelength \f rag ~ 47rif grows fastest, 



where H = c^/ttGE and c s = yksT '/ '/xmn is the sound velocity (Simon 1965). Although 
the primordial gas clouds we consider are not isothermal, the dispersion relations are 
weakly dependent on equation of state (see Fig.l of Larson 1985) and we use the isothermal 
dispersion relations as typical ones. As shown by Miyama et al. (1987a, b), an equilibrium 
sheet tends to fragment into filaments, whose line density is estimated by 

Since the above M is greater than M C (T) the filament formed by fragmentation of an 
isothermal sheet cannot be in an equilibrium configuration but begins to collapse. As 
shown by Inutsuka & Miyama (1992) such a cylinder with M > M C (T) does not fragment 
immediately. To estimate the epoch of the fragmentation of a collapsing cloud we 
should compare the dynamical time scale of collapse defined by td yn = p/-£ with that of 
fragmentation, tf rag . When td yn & tf rag , the cylinder collapses so fast that the density 
perturbation does not grow similar to the spherical case. The collapsing cloud fragments 
when the condition 

tdyn ~ tfrag: (8) 

is satisfied, i.e. when the collapse almost stops. In such a situation we can apply the 
perturbation analysis of the equilibrium sheet and the cylinder. 

In the case of the isothermal cylinder in equilibrium, the perturbation with wavelength 
Xfrag ~ 2irR grows fastest, where R is the scale radius of the cylinder defined by 
R = J Mj Trp with M and p being the line density and the central density, respectively. 
The e- folding time scale of the fragmentation is given by tf rag = 2.1/^2irGp , (Nagasawa 
1987). The mass of the fragment formed by the fragmentation of the isothermal cylinder is 
estimated by 

Mf ra g — Xf rag M ~ 2ttRM. (9) 



In our model we estimate the mass of the fragment by equation (j9|) at the epoch when the 
condition (Eft is satisfied. 
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4. The results of numerical calculations 

To obtain initial conditions of a cylindrical cloud which is formed by the fragmentation 
of the post shock cooled sheet, we recalculate a typical case of SK with z = 5, VL^h 2 = 0.1, 
and shock velocity v s = 200km s _1 . As initial conditions of our cylindrical cloud 
we adopt physical parameters of the post shock one dimensional flow at the point 
when the cooling time of the sheet (~ the dynamical time scale of the sheet) 
becomes equal to the time scale of the fragmentation, so that the sheet is expected 
to fragment into filaments. The initial condition of the cylinder thus obtained is 
n = 6.9cm" 3 , T = 330K, / Ha)0 = 2.3 x 10~ 3 , and / e -, = 2.7 x 10~ 4 . Under the above initial 
condition we calculated for the cases M/M C (T ) = 1.0, 1.3, 1.5, 1.7, and 2.0 (case A~E, 
see Table 1), for comparison. 

Since the equilibrium sheet tends to fragment into the cylinder of line density 
M ~ 2M C (T ) (eq. [I], |7]), we show the evolution for case E as representative, for which 
the mass of the fragment becomes minimum. The evolution on the n-T (number density - 
temperature) plane and virial temperature (eq. [[II] ) are shown in Figure [^a. From Figure 
[I]a we can see that the cylinder continues to collapse even after the cloud becomes optically 
thick to the line emission of hydrogen molecules (n ^ 10 10 cm~ 3 ). The fact that the cloud 
can collapse isothermally as a function of time by radiative cooling (n ^ 10 12 cm~ 3 ) is the 
characteristic feature of the cylindrical shape of the cloud, for which the gravitational force 
is proportional to Br 1 . Since temperature becomes constant, hydrogen molecules does not 
dissociate. Figure [l]b shows two time scales tf rag and t dyn normalized by the initial free-fall 
time scale (1/ \/2irGpo) and the mass of the fragment (eq. ||). Comparing two time scales, 
we conclude that the collapsing cylinder fragments at n ~ 3 x 10 12 cm~ 3 and the mass of the 
fragment becomes ~ 2Mq for case E. 

In Table [1] number density n, hydrogen molecule fraction /h 2 when the fragmentation 
occurs (column 3 and 4, respectively), and the mass of the fragment Mf rag (column 5) are 
listed for each line densities M (column 2). From Table [l] we can see the tendency that 
the larger the initial line density M, the smaller the mass of the fragment Mf rag and Mf rag 
converges to several solar mass for larger line densities. Interpretation of these results will 
be stated in the next section. 



5. Analytical estimate of the minimum fragment mass 

It is interesting to ask what determines the mass of the fragments of the collapsing 
cylinder for case E. The collapsing cylinder fragments when the condition (|8|) is satisfied. 
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As the cylinder collapses isothermally as a function of time finally, we can set de/dt ~ 
in equation (|6]) and from this we obtain t dyn ~ (7 — 1) t coo i, where 7 = 5/2 for diatomic 
molecular gas. Since the cloud cools by the line emissions of hydrogen molecules and is 
optically thick to these lines emission, t coo i can be estimated by 



knT 



tcool ~ Vrot^oc: (10) 

where a = 2ir 5 kj 3 /l5h 3 c 2 is the Stefan-Boltzmann constant, and a c is the effective number 
of line emissions. Since the optical depth is ~ 100 at most, the wing of line profile may not 
be important (Rybicki & Lightman 1979) and the line broadening is determined by Doppler 



broadening as Av/v = t>H 2 / c = ykpT /m^c 2 . The effects of chemical heating/cooling are 
not included in t coo i because the fraction of hydrogen molecule fu 2 is almost unity by the 
epoch of the fragmentation. The temperature of the cylinder is estimated by the virial 
temperature as 

k B T = ^m R GM. (11) 

From the above equations we obtain 

M frag ~ 2nRM ~ ,/__ « (12) 



where mp[ = y'hc/G is the Planck mass. Equation fl!2 ) means that Mf rag is essentially 
equal to the Chandrasekhar mass (~ mpjm^ here m p is proton mass). 

So far we did not consider other heating processes such as shocks and turbulence which 
likely to occur in the collapsing cloud. All these processes tend to halt the collapse, so that 
the epoch of fragmentation becomes earlier and, as a result, the mass of fragment increases. 
So equation ([12]) gives the lower limit on the mass of the fragment which is formed from the 
primordial gas cloud. Our estimate shows even the primordial gas cloud which cools most 
effectively by hydrogen molecules cannot fragment into a smaller mass than equation (|T2|). 



For smaller line densities the gas pressure force halts the collapse before three-body 
reactions (Palla et al. 1983) completely convert atomic hydrogens to molecules. In these 
cases equation (|12"D is not applicable because clouds are optically thin to line emissions 
for /h 2 ^ 10 _1 and chemical heating by H2 formation must be considered when hydrogen 
molecules are formed by three-body reactions. These clouds fragment at lower density, in 
other words at larger scale radius, than the density at which they reach the isothermally 
collapsing stage (corresponding to n ^ 10 12 cm~ 3 of case E). Since Mf rag is proportional to 
the scale radius at fragmentation, this leads to the larger fragment mass. 
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6. Discussions 

The opacity-limited hierarchical fragmentation scenario was argued by several authors 
(Low & Lynden-Bell 1976; Rees 1976; Silk 1977). According to Rees (1976), the minimum 
Jeans mass of the spherical cloud is given by 

Mf = ( kBT V /4 r 1 ' 2 1 m3p \ 

\m H c 2 J /i 9 / 4 ' 

where f (^ 1) is the efficiency of the radiation of the cloud compared with the black body 
radiation and it depends on the details of the cooling process and the opacity. As the 
present molecular clouds have various molecules and grains, f is nearly equal to unity when 
the cloud becomes optically thick to the radiation. In the case of the primordial gas cloud 
which is cooled by line emissions of hydrogen molecules, f is proportional to VT ( Doppler 
broadening ) and T-dependence of Mp completely vanishes. Apparently this gives the same 
conclusion as ours. However previous authors argued the spherical collapse, which does not 
lead to fragmentation, and our argument would be more relevant. 

Although our estimate is order of magnitude and the exact mass of the fragment cannot 
be estimated, the existence of the lower limit on the mass of the fragment we found would 
have significant meanings. For example, the following scenario for formation of our Galaxy 
may be possible: the existence of lower limit on the fragment mass cause the mass of the 
first stars to be massive, so that the first stars complete main-sequence stage within the age 
of the universe and evolve to dark remnants (e.g. white dwarfs, neutron stars, and black 
holes) in the Galactic halo after some mass loss, whose gas fall to the Galactic disk to form 
it. This scenario would be consistent with infall model (Lynden-Bell 1975; Vader & de Jong 
1981; Clayton 1988), which was proposed to explain G dwarf problem (the observed paucity 
of metal-poor stars, see Rocha-Pinto & Maciel 1996 for recent observation). Remnants 
thus formed by the above scenario may be MACHOs, whose estimated mass (~ 0.5Mq, 
Bennett et al.1996) is consistent with that of white dwarfs. It is very interesting that some 
of astrophysical problems, including non-existence of zero metal stars, may be explained in 
a consistent way by such a simple scenario. 

We would like to thank H. Sato for useful discussions, T. Chiba for critical reading, and 
S. A. Hayward for checking English. This work is supported in part by Research Fellowships 
of the Japan Society for the Promotion of Science for Young Scientists, No.6894(HU), 
3077(HS) and by Grant-in-Aid of Scientific Research of the Ministry of Education, Culture, 
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Table 1. initial line density and results of numerical calculations 



case 


M/M C (T ) 


n(cm 3 ) 


/h 2 


M frag /M Q 


A 


1.0 


1 x 10 8 


5 x 10- 3 


1 x 10 2 


B 


1.3 


1 x 10 9 


4 x 10~ 2 


4 x 10 1 


C 


1.5 


2 x 10 10 


4 x 10- 1 


2 x 10 1 


D 


1.7 


5 x 10 11 


9 x 10- 1 


3 x 10° 


E 


2.0 


3 x 10 12 


1 x 10° 


2 x 10° 







5 1( 
log 10 n (cm -3 ) 



Fig. 1. — a) The thermal evolution of a cylindrical cloud on the n-T plane (solid curve) for 
case E. The critical temperature is also plotted(dashed curve), b) Evolution of time scales 
and mass of a fragment for case E. Time scales are normalized by the initial free-fall time. At 
first the time scale for collapse, td yn (dashed curve), is shorter than that for fragmentation, 
tfrag (solid curve), by an order of magnitude. But finally tf rag becomes shorter than t dj/n , 
at which the collapsing cylinder is thought to fragment and the mass of a fragment (dotted 
curve) becomes ~ 2Mq. 



